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We present an analytical approach to the theory of nonlinear propagation in gases of femtosecond 
optical pulses with broad-band spectrum . The vector character of the nonlinear third-order polar¬ 
ization of the electrical field in air is investigated in details. A new polarization state is presented 
by using left-hand and right-hand circular components of the electrical field . The corresponding 
system of vector amplitude equations is derived in the rotating basis. We found that this system of 
nonlinear equations has 3D -I-1 vector soliton solutions with Lorentz shape. The solution presents a 
relatively stable propagation and rotation with GHz frequency of the vector of the electrical field in a 
plane orthogonal to the direction of propagation. The evolution of the intensity profile demonstrates 
a weak self-compression and a week spherical wave in the first milliseconds of propagation. 

PACS numbers: 42.65.Sf, 52.38.Hb 


INTRODUCTION 


When a femtosecond laser pulse with power above the 
critical for self-focusing propagates in air, a number of 
new physical effects are observed, such as long-range self- 


channelin; 
emission 
spectra 
ized pulse 



coherent and incoherent GHz and THz 
, asymmetric pulse shaping, super-broad 
3, polarization instability of linearly polar- 
polarization rotation [8|, self-compression 
0 and others. A remarkable effect is also that in lidar 
experiments the light filaments propagate over distances 
up to 9-10 kilometers in vertical direction, preserving 
their spectrum and shape 0. In a typical experiment 
in the near zone (up to 1 — 3 m from the source), when 
the initial pulse intensity exceeds I > lO^^IU/cm^, self- 
focusing and self-compressing start, which makes the kz 
spectrum broad band and asymmetric Ri kg. The 
process increases the core intensity up to lO^'^lT/cm^, 
where a short non-homogenous plasma column in the 
nonlinear focus is observed. Usually the standard model 
describing the propagation in the near zone is a scalar 
spatio-temporal nonlinear paraxial equation including in 
addition terms with plasma ionization, higher order Kerr 
terms, multiphoton ionization and others [13, El- The 
basic model works properly in the near zone because 
of the fact that the paraxial approximation is valid for 
pulses with narrow-band spectrum Akz << kg. At dis¬ 
tances longer than 10-20 meters from the laser source, 
where the stable filament is formed, plasma generation 
and higher-order Kerr terms are also included as nec¬ 
essary for the balance between the self-focussing and 
plasma defocussing and for obtaining long range self¬ 
channeling in gases. However, the above explanation of 
the filamentation is difficult to be applied at such dis¬ 
tances. As reviewed in lll-ll7j|. the plasma density at 
long distances from the source is too small to prevent 
self-focusing. There are basically three main character¬ 
istics which remain unchanged at these distances - the 


broad-band spectrum, the coherent GHz generation and 
the width of the core, while the intensity in a stable 
filament drops to a value of VdA^Wjcm}. The higher- 
order Kerr terms for pulses with intensities of the order 
of / ^ VS^’^Wjcrr? are also too small to prevent self¬ 
focussing. The experiments, where observation of long- 
range self-channeling without ionization was realized |12l - 
m , show the need to change the role of the plasma defo- 
cusing at such distances with another effect. In addition, 
there are difficulties with the physical interpretation of 
the coherent GHz radiation as a result of plasma gen¬ 
eration. The light pulse near the nonlinear focus emits 
incoherent and non-homogenous plasma, while the co¬ 
herent GHz radiation requires homogenous plasma with 
fixed electron density of the order of 10^® cm~^. Only ho¬ 
mogenous plasma can generate coherent GHz emission, 
but such kind of plasma is absent in the process of fila¬ 
mentation. In the real experiments with propagation of 
a single filament at distances more than 20 — 30 m from 
the source in air, the following basic characteristics are 
found: 

1. Broad-band spectrum [Akz ko). 

2. Intensity of the order of / ~ 10^^“^^ W/cm^. 

3. Absence of plasma at long distances. 

4. Asymmetric relatively stable (Lorentz) spectral and 
longitudinal shapes. 

5. Goherent GHz generation. 

Recently in [ 2 ^ we developed a scalar ionization-free 
non-paraxial nonlinear model, which gives the above 
characteristics of the stable filament. The analytical and 
the numerical results describe correctly the linear and 
nonlinear evolution of narrow-band and broad-band laser 
pulses. In addition it was found that the equation has 
exact Lorentz-type soliton solutions in approximation of 
neglecting the GHz oscillation. Still, this theory cannot 
resolve some difficulties. The main problems are: 

1. Peak instability of the soliton solution under small 
initial perturbations. 
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2. The soliton solution is obtained after neglecting the 
GHz oscillation. 

3. The soliton solution has one free parameter. 

4. There are problems with the conservation law of the 
nonlinear operator when we use the GHz oscillation. 

To solve the above problems in this paper, we propose 
a nonlinear vector generalization to the model. 


NONLINEAR POLARIZATION 

The self-action process broadens the pulse spectrum - 
starting from a narrow-band pulse, the stable filament 
becomes broad-band far from the source. In recent pa¬ 
pers [13, it is shown that the evolution of broad-band 
pulses like filaments can be described correctly by useing 
the generalized nonlinear polarization 


P"' = n2 E, (1) 

which includes additional processes associated with third 
harmonic gen eration (THG). A more precise analysis pre¬ 
sented in [2l| demonstrates that the polarization of the 
kind (ID) is not applicable to a scalar model, because the 
corresponding Manley-Rowe (MR) conservation laws are 
not satisfied. That is why we investigate two-component 
electrical vector field. The generalized nonlinear polar¬ 
ization m is quite simple in terms of left-hand and right- 
hand circular components. Let us now present the elec¬ 
tric vector field of a pulse as a linear decomposition of 
left- and right-hand circular complex components. 


AE± - 


1 d^E+ Att 


dt"^ dE 


I J R^'^\T)E±{t — T)dT 


(7) 


poo poo poo 

+ / / R^^\T,T,T)[E±{t-T)f E^{t-T)dT^\ 

Jo Jo Jo ■’ 


To obtain amplitude equations, we use in mind also the 
causality principles (no negative time r ) to the response 
functions and their Fourier presentations (w) and 


poo 

= / R^^\t) exp{iujT)dT 

Jo 


( 8 ) 


poo poo poo 

X^^\oj,uj,uj) = / / / i?^^^(T, r, r)exp(3fwr)(ir^. 

Jo Jo Jo 


Reduction of the integrals ([5]) from 0 to infinity is equal 
to cosFourier transforms with the properties 


(9) 

X^^\uj,uj,uj) = x''^\-0J,uj,uj) = x^^^(-u;,-a;,a;)... 

The key point in the following transformations is the 
fact that the circular components of the electrical field, 
as well as the amplitude functions, are orthogonal in the 
complex plane. That is why their Fourier presentations 
are also written in orthogonal basis 


E{x,y,z,t) = E+{x,y,z,t)d+ + E_{x,y, z,t)d- 

(2) 

E+{t -'^) = ^ 

pOO 

/ Ej^[(jj)exp[—i(jj(t — T)du} 

(10) 

where the circular- polarization unit vectors are 



— OO 


a± = (x± iy)l'/2. 

(3) 

1 

1 

II 

poo 

/ E-{(jj) exp{iuj{t — T)duj. 

J —OO 

(11) 


If we now represent P"* in terms of its circular compo¬ 
nents as 


P = P+d+ + P_tT_, 

we find that the components are given by 

P+ = 2n2 (P+P-) 

P_ = 2n2 {E‘iE+) . 


(4) 


(5) 

( 6 ) 


Substituting (Uni) - (ED into the right-hand side and also 
to the last term of left-hand side in dD, and using the 
spectral properties ([9]) of the response functions (|8]), the 
wave system can be written as 


AP+ = - 


k'^{uj)E+{uj) ex'p{—iLjjt)dijj 


( 12 ) 


pOO pOO pOO 

' —oo J —oo J —oo 




E, 


E- exp{—iojt)du)^ 


BASIC SYSTEM OF EQUATIONS 

The decomposition ([D allows us to rewrite the non¬ 
linear vector wave equation in the following system of 
equations. 


AP_ = - 


k‘^{uj)E-{uj) exp{iujt)duj 




exp{iiot)du!^, 


( 13 ) 
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where k‘^{oj) = w^e(a;)/c^; = 

k'^{uj)n 2 {uj) are the square of the linear and nonlinear 
wave vectors and 


£(a;) = 1 + 47rx^^^(a;) (14) 

n2(w) = 47rx^^^(a;,a;,a;)/e(a;). (15) 

Lets us now introduce complex amplitude functions using 
the substitutions 


E+{r,t) = A+{r,t) exp[i {koz - ojot)] (16) 


E- (r, t) = A- (r, t) exp [—* {k^z — wq!)] . (17) 

where r = {x,y,z), ojq is the carrier frequency and ko 
is the linear part of the wavevector ar the carrier fre¬ 
quency. Applying the translation theorem to the Fourier 
presentations of the electrical field components we have 


E^(r,uj) = exp(ikoz)A+(r,uj — uJo) (18) 

E-{r,uj) = exp(—iA:oz)A_(r, w — wq). (19) 
Substituting (IT6l) - (fT^ into Eqs. m we obtain 


A . .A., dA+ , , , 

Ay4_|_ + 2 z/po—-^0-^+ “ 

az 

k'^{uj)A+ exp [—i{w — wo)t] t)du} 


( 20 ) 


°° r A -| 2 .. 

k^i{uj) A+ A-exp[—i{uj — uJo)t]duj^ 


dA- 

AA_ — 2ikQ— -= 

oz 

k‘^{uj)A- exp [i{(jj — wo)^] t)duj 


( 21 ) 


fcn/(w) 


A_ 


A_|_ exp [i{oj — Wo)!] du)^ 


We note that all functions in the integrals on the right 
side of Eqs. (l20ll and ([ST]) depend form on the frequency 
difference (w — wq), except for the linear k‘^{iu) and non¬ 
linear A:^;(w) wave vectors. Then we expand k'^{uj) and 
k'^i(uj) as a power series of the same difference (w — wq) 


fc^(w) = fcp + 2—^ (w — Wo) 

Vgr 

+ (~ 2 —^ kQk''\ (w — Wo) -|- ... (22) 

\^gr J 


= fc^(w)n2(w) = fco^2U=CA 


2ko A , ,2^’^2| 

-n2(wo) -l- kQ— — \u]=ujo 

Var OOJ 


(w — Wo) -I- ..., (23) 


where Vgr is the group velocity, and k" is the group ve¬ 
locity dispersion. After replacing /c^(w) and k'^i{uj) with 
their series (l2^ - (l23t . all functions in the integrals from 
the right side of Eqs. o and m depend on the fre¬ 
quency difference and we can integrate over all values of 
(w — Wo). Thus, we obtain the following system of ampli¬ 
tude equations in circular basis (up to second order of the 
group velocity dispersion and zero order of the nonlinear 
dispersion). 


2iko 




1 dA+\ 

Vgr Ot J 


1 

Vgr dt J 


= AAi - 


1 -f p d^A^ 
dt"^ 


gr 

r,2 


-|-2fcQn2A^A_ 


(24) 


= AA_ - 


1 -H /3 d'^A_ 
v'jr 

4-2/c Q u 2 A _ A _|_, 


where A is 3D - (x, y, z) Laplace operator, and /3 is a 
number, connected with the dispersion characteristics of 
the medium (/3 = koVgrk”). We note here that in gases 
the dispersion is weak and the series (HH) are strongly 
converged up to single cycle regime (broad-band pulses). 
Thus, the non-paraxial system of equations (l24l) describe 
correctly the evolution of laser pulses in gases up to a 
single cycle regime. It is important to mentoin that from 
Eqs. (1241) paraxial spatio-temporal ap pro ximation can be 
derived for narrow-band laser pulses 22| only. The fila- 
mentation experiments demonstrate quite different pulse 
evolution: the initial laser pulse (!o > 50/s) possesses a 
relatively narrow-band spectrum { Ak^ <C A:o) and during 
the process the initial self-focusing and self-compression 
the spectrum broadens significantly. The broad-band 
spectrum {Akz ^ kg) is one of the basic characteristics 
of the stable filament. That why we do not reduce more 
Eqs. and try to solve them for the case when the 

pulse has a large spectrum. 

Another standard restriction in the filamentation the¬ 
ory is the use of one-component scalar approximation of 
the electrical field E. This approximation, though, is 
in contradiction with recent experimental result^ where 
rotation of the polarization vector is observed [8|. For 
this reason in the present paper we use non-paraxial vec¬ 
tor model in circular basis (I24L in which the nonlinear 
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effects are described by the nonlinear polarization com¬ 
ponents The dispersion number in air is very small: 
P = kQVgj.k" ~ 2.1 X 10“^, so we can solve Eqs. (1^ in 
approximation up to the first order of dispersion. Addi¬ 
tionally, we will use the normalized amplitude functions 
A± = AoA± to rewrite the system (0^ in the form 


2iko 



—2iko 



1 

vj, dt^ 

(25) 

1 d^A_ 

Vgr dt J df^ 

+lA^-A+, 


dt 


1 dA_ 


— A A 


where 7 = 2 n 2 A:gAQ is a nonlinear coefficient. 


VECTOR SOLUTION AND VECTOR ROTATION. 
THE FILAMENT AS A WEAK ROGUE WAVE 

The nonlinear system of equations ((25ll has exact soli¬ 
tary vector solution when 7 = 2 and the spectral width 
of the pulses reaches the value Ak^ « fco 


E^{x, y, z, t) = 

(29) 

2 

Ey{x, y, Z, t) = cos [Afc^ {Vph - Vgr) t]. 

The 3D -I- 1 Lorentz type solution (ESI), presented in 
Cartesian coordinates, gives oscillation of the electrical 
vector E = {Ej;,Ey,0) in the {x,y) plane. It can be 
seen directly that the frequency of oscillation is equal to 
the carrier to envelope frequency cocef = ^0 (vph — Vgr). 
The corresponding longitudinal spatial carrier to enve¬ 
lope wave number is kcef = uJcef(vgr- Detailed investi¬ 
gation on the evolution of the solution (l29l) at distances 
more than 10 — 20 meters from the initial point, reveals 
weak self-compression. That is why we consider the fila¬ 
ment in this zone as closer to a weak Rogue wave. The 
evolution of the profiles of the electrical field components 
3?(ifa;)-Fig. la and K(ify)-Fig. lb. are plotted in Fig. 1. 
The periodical exchange of energy between components, 
due to nonlinear mechanism, leads to rotation of the elec¬ 
trical field vector in a plane orthogonal to the direction 
of propagation, with time period Tcef ~ 1 — 2 x 10“^° s 
and spatial period Ace/ — 3 — 6 cm. 

CONCLUSIONS 


A+{x, y, z, t) = exp [f A/c^ (z - VgA)] 

(26) 

2 

A_ (x, y,z,t) = -r—-ZT exp [-iAk^ (z - Vgrt )], 


where f = ^x'^ + y'^ + {z - ikcefY - - ikcef/vgrY 

and kcef is determinated below. The solution of the 
corresponding vector electrical field can be written af¬ 
ter multiplying the amplitude functions (j26l) by the main 
phases dTi-dni) 


E+ {x, y, z, t) = exp [iAk^ {vph - Vgr) t] 

(27) 

2 

A. (a;, y, z, t) = ^^ 7^2 “ '^ 9 ^) P ■ 

Let us turn from the left-hand and right-hand circular 
components ([3]) to the standard Cartesian coordinates 


The starting point of our investigation in this paper 
is the fact, that the generalized nonlinear polarization 
o arise new polarization state in circular basis. This 
polarization state leads to periodical exchange of energy 
between the electrical components E^ and Ey of a laser 
pulse. To derive the corresponding amplitude equations 
associated with this polarization, we take into account 
the orthogonality in the complex plane of the left-hand 
and right-hand circular components. We find that the 
obtained system of amplitude equations (1251) has exact 
(3D-I-1) Lorentz type soliton solutions (l26l) . Our soli- 
ton solution is obtained for pulses which satisfy the addi¬ 
tional condition Akz ~ fcg. The diffraction of broad-band 
pulses is not the Fresnel one 23|, [2^ , which leads to the 
conclusion that the soliton appears as a balance between 
semi-spherical (Fraunhofer type) diffraction and nonlin¬ 
ear self-focusing. The solution gives also a rotation of the 
vector of the electrical field with the carrier to envelope 
frequency. 


E, = {E+ + E_)lV2, Ey = {E+ - E-)/{iV2). (28) 

The solution (l?7l) written in Cartesian coordinates has 
the form 
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FIG. 1. Evolution of the profiles of the electrical field components ?R.{Ex)-Fig. la and K(i?y)-Fig. lb of the solution (1291) . The 
periodical exchange of energy between components, due to nonlinear mechanism, leads to rotation of the electrical field vector 
in a plane orthogonal to the direction of propagation. 
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